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A simple quantum mechanical model consisting of a discrete level resonantly coupled to a con- 
tinuum of finite width, where the coupling can be varied from perturbative to strong, is considered. 
The particle is initially localized at the discrete level, and the time dependence of the amplitude 
to find the particle at the discrete level is calculated without resorting to perturbation theory and 
using only elementary methods. The deviations from the exponential decay law, predicted by the 
Fermi's Golden Rule, are discussed. 
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The transition of a quantum particle from an initial 
discrete state of energy e into continuum of final states 
is considered in any textbook on quantum mechanics. It 
is well known that perturbation theory approach, when 
used to solve the problem, leads to Fermi's Golden Rule 
(FGR), which predicts the exponential decrease of the 
probability to find the particle in the discrete state. It 
is also well known that, even for a weak coupling be- 
tween the discreet state and the continuum, this result 
(exponential decrease of probability) has a finite range of 
applicability, and is not valid either for very small or for 
very large time (see e.g. Cohen- Tannoudji et. al. [l[). 
This complies with the theorem proved 50 years ago, and 
stating that for quantum system whose energy is bounded 
from below, i.e., (0, oo) the exponential decay law cannot 
hold in the full time interval 0, S 0] ■ The same state- 
ment remains valid when the discrete state is coupled to 
the continuum bounded both from below and above, the 
model we presently consider. Fortunately, the model is 
exactly solvable; the solution in the frequency represen- 
tation obtained within the Green's functions (GF) for- 
malism is presented in where the model is called the 
Fano- Anderson model). In this paper we solve the prob- 
lem of tunneling in the framework of Quantum mechan- 
ics proper We concentrate on the calculation of the time 
dependent non-decay amplitude, the stage which is typ- 
ically not given proper attention to within the GF for- 
malism [1]. We analyze the relation between the exact 
results and those given by the FGR. 

The system consists of the continuum band, the states 
bearing index k, and the discrete state d. The Hamilto- 
nian of the problem is 

if = ^ c^fe Ifc) (fcK e M) (d| + ^ (Ffc |fc) {d\ + h.c) , (1) 

k k 

where \k) is the band state and \d) is the state localized 
at site d; h.c. stands for Hermitian conjugate. The wave- 
function can be presented as 



m=9{t)\d)+Y.^{k,t)\k), 



(2) 



Schroedinger Equation for the model considered takes 
the form 



^-M = eg{t) + Y^V:h{k,t) 



dt 

db{k,t) 

I ; 

dt 



= LJkb{k,t) + Vkg{t) 



Making Fourier transformation (Im uj > 0) 



g{t)e'^'dt, 



(3) 



(4) 



we obtain 



-i + ujg{uj) = eg{uj) + '^Vkb{k,uj) 

k 

ujb{k,uj) = LLikb{k,uj) + Vkg{oj)- 



(5) 



For the amplitude to find electron at the discrete level, 
straightforward algebra gives 



2ni J 



where 



s(^) 



1 



OJ - 



E 



e - S(w) 
\Vk? ^ 

UJ -LUk' 



(6) 

(7) 

(8) 



and integration in Eq. ([6]) is along any infinite straight 
line parallel to real axis in the upper half plane of the 
complex to plane. Notice that if, in addition, we define 
g{t) = for t < 0, then g{t) is just the (single particle in 
an empty band) GF in time representation, and g^io) is 
the same GF in frequency representation. The quantity 
T,{u!) in the GF formalism is called self-energy (or mass 
operator). 

For tunneling into continuum, the sum in Eq. ([S]) 
should be considered as an integral, and Eq. ([8]) takes 
the form 



with the initial conditions a(0) — 1, b(k,0) = 0. 



Et 



A{E)_ 



dE, 



(9) 
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FIG. 1: Contour used to evaluate integral (|6}. Radius of the 
arc goes to infinity. 



where 



(10) 



where and the hmit of integration are the band bottom 
Ell and the top of the band Et . We would like to calculate 
integral ^ closing the integration contour by a semi- 
circle of an infinite radius in the lower half-plane. Thus 
we need to continue analytically the function g{uj) which 
was defined initially in the upper half plane (excluding 
real axis) to the whole complex plane. We can do it 
quit simply, by considering Eqs. ([7]) and ^ as defining 
propagator in the whole complex plane, save an interval 
of real axis between the points Ei, and Et , where Eq. ([9]) 
is undetermined. (Propagator analytically continued in 
such a way we'll call the standard propagator.) Thus the 
integral is determined by the integral of the sides of the 
branch cut between the points E^, and E^. The real part 
of the self-energy T,' is continuous across the cut, and the 
imaginary part S" changes sign 



T."{E + is) = j:"{e- 



TrA{E) 



So the integral along the branch cut is 



A(i;)e- 



E, [E-e-i:'{E)Y +^^/^-'{E) 



Thus we have 

and the survival probability p{t) is 
P{t)^\9{t)\\ 



-0. (11) 



(12) 



(13) 



(14) 



In the perturbative regime |E'(e)|, |S"(e)| ^ e — 
Ef,, Et—e the main contribution to the integral IT^ comes 



from the region E 
sented as 



e. Hence the integral can be pre- 

A{e)e-'^'dE 



(15) 



and easily calculated to give the well known Fermi's 
golden rule (PGR) 



Pit) 



where 



l/r = 27rA(e) 



(16) 



(17) 



However, even in perturbative regime, the FGR has a 
limited time-domain of applicability [l| . For large t the 
survival probability is determined by the contribution to 
the integral ((T^) coming from the end points. This con- 
tribution can be evaluated even without assuming that 
the coupling is perturbative. Let near the band bottom 
(the contribution from the other end point is similar) 
A(£;) {E~ Eb)/^, where /? > 0. Then for large t 



-(/3+1) 



(18) 



The similar contribution comes from the top of the band. 

If there can exist poles of locator ([T]), we should add 
to the integral (fT2l) the residues 



9it)=I,utit) + Y.Rj, 
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(19) 



where the index j enumerates all the real poles Ej of the 
integrand, and 



iEit 



Rj - 



1 



dE I £=£;,■ 



(20) 



Notice, that the poles correspond to the energies of bound 
states which can possibly occur for E < Et or E > Et, 
and which are given by the Equation 



E,=e^ 



E 



UJ - UJk 



(21) 



If we take into account that normalized bound states are 

- Vk 



\E, 



i + E 



k (eP^ 



1/2 ' 



(22) 



then the residue can be easily interpreted as the ampli- 
tude of the bound state in the initial state times the 
evolution operator of the bound state times the ampli- 
tude of the state \d) in the bound state 



R, = {d \E.i) {E, \d) e 



iEit 



(23) 



3 



If the locator has one real pole at Ei, from Eq. (|T9)) 
we see that the survival probability p{t) when 
t — + cxD. If there are several poles, this equation gives 
Rabi oscillations. 

Notice that Eq. is just the well known result Q 



where 



-Im [g{E + is)] 



(24) 



(25) 



is the spectral density function. The first term in Eq. 
(I19p is the contribution from the continuous spectrum, 
and the second term is the contribution from the discrete 
states. 

Before proceeding further, consider two simple models. 
First model is a site coupled to a semi-infinite lattice. 
The system is described by the tight-banding Hamilto- 
nian 

H = -^Y.{\n) (n+l| + |n + l) (n\) 



-e\d) {d\-V{\d){l\ + \l){d\), 



(26) 



where \n) is the state localized at the n-th site of the 
lattice. The band (lattice) states are described by the 
Hamiltonian 



Hq = — cos k\k) {k\ , 




FIG. 2: Survival probability as a function of time for Ao = 
.02. 




(27) 
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FIG. 3: Survival probability as a function of time for Ao — .1. 



where |fc) — \/2 X^n sin(fcn) \n). Hence we regain Hamil- 
tonian ([T]) with Vfc = —\/2V sink; also 



^ r Ao(£;-sign(£;)/B^) 



|£;|>i 



0, 



\E\ > 1 



Ao\/r~B^, 1^1 < 1 ' 



(28) 



where Ao = 2V^. If we consider the case e = and 
I^^/2 < I the locator does not have real poles, so Eq. 
(fT2| after substitution of the results of Eq. (|28p gives the 
amplitude of the non-decay amplitude [6|]: 



dQ exp(ji cos Q) 



I - exp{~2iQ) 
1 + exp(-2i(5) 



,(29) 



where = 1 — Aq. 

For numerical calculations we consider the second 
model: 



A{E) = Ao = const. 



(30) 



The units and zero of energy we'll chose such that Ef, 
— l,Et — 1 ; thus we get 



J:'{E) = Aoln 



E + 1 



E- 1 



(31) 



There are two real poles of the locator, given by the Equa- 
tion 



E -e - Aoln 



E+1 



= 0. 



E- 1 

The time we'll measure in units of the FGR time t 
1/t ^ 2nAo. 



(32) 



(33) 



For the sake of definiteness we'll chose e — —A. For 
Aq = .02 (see Fig. 2) we observe the FGR regime, say, 
up to t = 9. For A = .1 (see Fig. 3) the FGR regime 
is seen up to t = 3. For Aq = .2 (see Fig. 4) the FGR 
regime is absent. The Rabi oscillations we see already 
at Fig. 3 and still more vividly at Fig. 4. The FGR is 
not valid for small t either. (From Eq. ^ it is obvious 
that the expansion of g{t) is g{t) = 1-1- kt^ -t- . . . , which 
gives quadratic decrease of the non-decay probability at 
small t.) The region of small t, shown at Fig. 5, clearly 
demonstrates quadratic decrease of p with time. 

The standard analytical continuation for the propaga- 
tor, which corresponds to using Eq. ([9|) in the whole u> 
complex plain (save a part of real axis, where this equa- 
tion is not determined) gives only one sheet of g, which is 
in fact a multivalued function of lu. An indication to the 
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FIG. 4: Survival probability as a function of time for Ao = .2. 
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FIG. 6: Alternative way to analytically continue the propa- 
gator into the lower half-plane. The arrows show the curves 
of analytical continuation in between the cuts and outside. 
Radius of the big arc and length of the cuts go to infinity. 
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FIG. 5: Survival probability for Ao 
process of tunneling). 



.02 (beginning of the 



existence of other sheets we obtain from further analysis 
of what seems to be the (almost) trivial result: FGR at 
perturbative regime. In fact, in this regime, Eq. (|16p we 
could obtain directly from Eq. ([6]), changing exact Green 
function ([7]) to an approximate one, which may be called 
the FGR locator 



1 



tj-e-S'(e)+i7rA(e) 



(34) 



(Notice, that whichever approximation we use for S(ti-'), 
the property a{t = 0) = 1 is protected, provided E does 
not have singularities in the upper half-plane.) Thus ap- 
proximated, locator has a simple pole oj = e -I- S(e), and 
the residue gives Eq. (fTB|) . However, the locator which is 
used to obtain Eq. (p5)) and the FGR locator have totally 
different singularities. So the fact that the locators give 
the same survival probability (even in finite time inter- 
val)demands explanation, which is as follows. The points 
Eb and Et are the S(cl') branch points. Hence the values 
of the self-energy in the lower half-plane depend upon 
the curve along which we continue the function from the 
upper half-plane 0]. The standard result corresponds to 
continuing Ti(uj) along the curves which circumvent the 
right branch point clockwise and the and the left branch 
point anti-clockwise. On the other hand, we could use 



a different continuation, along the curve passing through 
the part of real axis between Ei, and Et. In this case, 
Eq. (j34p which is valid in the perturbative regime near 
the point a; = e in the upper half-plane is valid in lower 
half-plain also, thus giving a pole in a different sheet of 
the multivalued propagator. More specifically, to obtain 
a single sheet, which we need to calculate the integral ^ 
we could have made the cuts from the branch points to in- 
finity and continuing the function between the cuts along 
the curves passing through the part of real axis between 
Eh and and outside as we did it previously. This way 
to make analytic continuation, and hence to calculate the 
integral ^ is presented on Fig. 6. ( The treatment of 
tunneling from the discrete level into a semi-bound con- 
tinuum presented in the paper by Onley and Kumar 
corresponds, in fact, to an analytic continuation similar 
in spirit, to that presented on Fig. 6.) Notice, that be- 
cause of the exponential decrease as a function of t of the 
integrand in the cut integrals appearing in this analytic 
continuation, in contrast to oscillatory behavior of the 
real axis cut integral, such analytic continuation is more 
convenient for the numerical calculations of the large t 
behavior of the non-decay amplitude. 

For example, for the model with constant A(_E), thus 
continued self energy is 



S(w) = Ao log 



1 



1 



(35) 



where log is defined as having the cuts [— 1,— 1 — ^cxd) and 
[1,1 — zoo), and the phase — tt at the real axis between 
— 1 and 1. For the model of the semi-infinite chain the 
relevant self energy sheet is 



E(w) = Ao (t^ - - l) 



(36) 



where the square root is defined as having the phase 7r/2 
at the real axis between —1 and 1. In both cases the 
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result for g{t) would include the integrals along the cuts 
presented at Fig. 6. 

A very interesting problem would be the analysis of 
the whole Rimann surface of the propagator, but this is 
unfortunately, at least at the present moment, above the 
author's head. 

CONCLUSIONS 

In this paper we solve the problem of tunneling from 
a discrete level into continuum. We show, how the ba- 
sic notion of the GF formalism, like frequency dependent 
discrete level propagator, self energy and spectral density, 
appear within the basic quantum mechanics. We concen- 
trate on the calculation of the time dependent non-decay 
amplitude, the stage which is typically not given proper 
attention to within the GF formalism [5|. We show that 
the exponential time dependence of the non-decay prob- 
ability given by the FGR is an approximation valid in 
perturbative regime and only for intermediate times. The 
large time dependence of the non-decay probability de- 
pends crucially upon the details of the band structure and 
hybridization interaction. We also look closely at those 
analytic properties of the propagator in the complex uj 
plane, which often pass unnoticed. 
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